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Abstract 

Maldacena has shown that the wavefunction of the universe in de Sitter space 
can be viewed as the partition function of a conformal field theory. In this paper, 
we investigate this approach to the dS/CFT correspondence in further detail. We 
emphasize that massive bulk fields are dual to two primary operators on the boundary, 
which encode information about the two independent behaviors of bulk expectation 
values at late times. An operator statement of the duality is given, and it is shown 
that the resulting boundary correlators can be interpreted as transition amplitudes 
from the Bunch-Davies vacuum to an excited state in the infinite future. We also 
explain how these scattering amplitudes can be used to compute late-time Bunch- 
Davies expectation values, and comment on the effects of anomalies in the dual CFT 
on such expectation values. 
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1 Introduction 

In p], it was proposed that the Hartle-Hawking wavefunction [2j in de Sitter space is 
equal to the partition function of a conformal field theory on the boundary. If such 
a dual field theory exists, then it gives a powerful method for computing in-in expec- 
tation values of bulk fields in the Euclidean vacuum at late times: instead of doing a 
path integral over the standard Schwinger-Keldysh contour, one can directly evaluate 
the wavefunction by computing the CFT correlators. Expectation values can then be 
calculated by squaring the wavefunction, inserting a product of bulk fields, and inte- 
grating over the CFT sources. 

This procedure was investigated in [I] in the context of a model of inflation with 
only massless fields. The massless scalar field operator in de Sitter space has only one 
independent behavior near future null infinity X + that can be identified with a primary 
operator in a CFT, implying that a massless scalar in the bulk corresponds to only 
one primary operator on the boundarj0. In this paper, we will generalize the results of 
[I] to the case where the bulk field is massive, and has two leading behaviors that can 
be identified with primary operators. This leads to a CFT correlator structure that is 
richer than in the massless case, since bulk correlation functions have more than one 
independent behavior near I + . 

Several important aspects of the duality were recently clarified in [4j. In particular, 
it was shown that extrapolating bulk correlation functions to the boundary does not 
give the CFT correlation functions defined in [JJ. One of the goals of this work is 



In fact, one can identify both independent behaviors of a massless bulk scalar field with CFT operators, 
but one of the operators has logarithmic correlators and is therefore not primary [3]- 
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to represent the correlators of [T] as extrapolated transition amplitudes between the 
Bunch-Davies vacuum and a certain final state, with a time-ordered product of fields 
inserted in between. This allows for an operator formulation of the correspondence, 
where the boundary operators are defined by pulling the bulk operators to the bound- 
ary, as in AdS/CFT [5]. 

The rest of this paper is organized as follows. In Section [21 we define the CFT oper- 
ators in terms of the boundary behavior of the bulk fields, and express the correlators 
of pQ as the boundary values of transition amplitudes. We also prove that modifying 
the final state in the amplitudes is accomplished by deforming the CFT by a certain 
marginal operator. In Section O we reinterpret the correlation functions as a scatter- 
ing matrix, and derive a simple diagrammatic expansion for this S-matrix. We make 
contact with the in-in formalism in Section U and show that Bunch-Davies expectation 
values at late times are computed by integrating over two sources in the dual CFT. 
Section [5] then concludes this work by discussing the effects of the conformal anomaly 
in the dual field theory on bulk expectation values. We argue that such anomalies are 
only visible in expectation values involving time derivatives of the bulk field <j>. 

2 Correlation Functions 

The goal of this section is to express the correlation functions of pQ in terms of bulk 
transition amplitudes. For simplicity, we will consider a scalar field (j) in the bulk, 
governed by the action 




(2.1) 



We work in the inflationary patch of dS^+i, where the metric reads 



ds 2 




OO < 7] < 0. 




Solutions to the equations of motion behave near X + as |6j 



4>(V,X) ~ 4> + {x)(-7]) A + 4>^(x)(-7]) 



(2.3) 



where the dimensions are 



We restrict our analysis to the case where A_|- are real, which is equivalent to only 
considering sufficiently light scalars. 

In anti-de Sitter space, the free field operator is expressed as a sum over normalizable 
modes, which are fast-falling at the boundary (see [7], for example). On the other hand, 
the field operator in de Sitter space has two behaviors at X + , and each of these can 
be identified with a CFT operator on the boundary. To be precise, we postulate the 
operator relation 



Here we have chosen a convenient normalization for 0±, which agrees with the AdS/CFT 
normalization up to a factor of ±i [5j[8]. Note that 0+ and 0_ are anti-Hermitian. 

It is important to emphasize that we are not assuming the existence of a dual 
CFT. Rather, the claim is simply that expectation values of 0± between bulk vacuum 
states obey the same constraints as correlation functions in a CFT. Also, note that 
our postulate (2.5) differs in spirit from much of the recent literature, which advocates 
a dS/CFT dictionary that is the analytic continuation from AdS/CFT [HH1[9]. One 
could imagine fixing future boundary conditions on a bulk field in de Sitter space by 
analytically continuing the field operator from AdS, but we instead find it convenient 
to fix boundary conditions by fixing the final state. The relation between the two ap- 
proaches will be clarified below; it turns out that the correlation functions of [U HJ E] 
are equal to the expectation values of + between two certain bulk vacua. The oper- 
ator O- acts as the source in the CFT of [HIUE], and captures the leading behaviors 
of expectation values. 

By dimensional analysis, we expect that 0± has conformal dimension A±. To 
confirm this, let us expand the operator <3? in terms of Bunch-Davies creation and 
annihilation operators |13j . 



$(77, x) 



±-0-(x)(- v ) A - + ±-o+(x)(-t,) 



(2.5) 



$(77, A?) 



V^{-ri) d/2 




(2.6) 
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where H u is a Hankel function of the first kind. Taking r] — > 0, we find 



0-{k) 



o+{k) 




(a) {-k) - a(k)) 



a\-k) -e iuw a(k)) 



(2.7) 



(2.8) 
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and these indeed transform as operators of dimension A±. In deriving the expressions 
(|2.7p and (|2.8p . we have assumed that v is neither an integer nor a half- integer - these 
cases, which include the case of a massless field, must be treated separately. 

Now let us define a state |+) that is annihilated by the slow- falling modes near 
J+, so that = [ID]. Suppose that we are interested in computing transition 

amplitudes from the Bunch-Davies vacuum |BD) to (+|, with a time-ordered product 
of $ fields inserted. These quantities can be computed via a path integral over fields 
obeying fixed boundary conditions in the infinite past and future, 

/ -<«r,=0M-r ? ) A + 

<+|T*(xi)---$(z„)|BD)= / [^MxO-.-^exp^]), (2-9) 

where we have absorbed an overall normalization into the path integral measure. The 
boundary conditions at early times require that <fi is purely negative frequency with 
respect to |BD) [I], while the late-time boundary conditions ensure that <fi is fast-falling 
near the boundary. 

In [H H] , it was shown that the infrared wavefunction in AdS analytically continues 
to the Bunch-Davies wavefunction in de Sitter space. Similarly, we see that the matrix 
elements (|2.9p are the analytic continuations of bulk correlation functions from AdS 
(related claims have appeared in [HE]). Indeed, the boundary conditions in (|2.9p 
analytically continue to the AdS boundary conditions under r/ = iz. Another way to 
see this is by computing the Green's function (+|T<3?(x)<3?(y)|BD), since the amplitudes 
(|2.9p are evaluated in perturbation theory using this Green's function. Using (|2.6p and 
the fact that C_|+) = 0, we find 

<+|T$(t7, £)$(»/, fc')|BD) = ~( m ') d/2 [j u (rjk)H v {rfk)9(f]-rf) 

+H u ( V k)J u ( V 'k)6( V ' - 77)] (2ir) d 5 d (k + k'), (2.10) 

which is indeed the analytic continuation of the bulk-to-bulk propagator from AdS [14J. 
This explains why the AdS propagator does not analytically continue to the two-point 
function in any de Sitter invariant vacuum, as was shown indirectly in [15^ 116] . The 
quantities (j2.9f) were referred to in [IT] as correlators with future boundary conditions. 

By taking the arguments of the amplitudes (|2.9p to the boundary and using the 
operator relation (|2.5p . we may compute the matrix elements of 0± between (+| and 
|BD). In fact, we claim that the matrix elements of 0+ are precisely the CFT corre- 
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lators defined in pQ. To see this, let us define the state 

10+) = exp (- J £h- d ^ + (k)0 + (-k)^ |+>, (2.11) 

where (ft+(x) is an arbitrary real function. Using the commutation relation 

[0+{k),0 -{%')} = 2iv(2ir) d 6 d (k + k'), (2.12) 
one may check that \4>+) is an eigenstate for 0_, 

0_(fc)|0+> = 2u/0+(fc)|0 + ). (2.13) 
The generating functional for matrix elements of C± between (+| and |BD) is therefore 

Z C ft[0±] = (+|exp yj ^L^O+i-k)^ exp j |BD) 

= exp hi* y j * BD [(-r?c) A -0+], (2.14) 



where we used the fact that $(r/, x) ~ — ^(— ry) A_ C_(x) near X + . Here we have 
introduced a late-time cutoff at 77 = rj c - When the source 0_ is turned off, we find 
that the argument of the Bunch-Davies wavefunction acts as a source for an operator 
of dimension A+, as shown in pp. 

Let us now comment on the effect of changing the final state (+| to a vacuum state 
that is instead annihilated by a linear combination of C + and O-. Up to an overall 
normalization, these vacua can be expressed as the squeezed states [15] 



|c> = exp (^c j ^ a \k)a\-k)\ |+>, (2.15) 

where c is a constant. Using (|2.7|) and (|2.8|) . we see that taking |c) instead of |+) as 
the final state is equivalent to inserting the exponential of a linear combination of the 
terms 

J d d xO + (x)0^(x), J d d xd d y^0^0- (2.16) 

into the generating functional (|2.14p . In fact, the operators that are integrated over x in 
(|2,16p are all marginal, so the resulting conformal field theory is a marginal deformation 
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of the original CFT, as conjectured in [15jo Note, however, that the deformation is 
nonlocal, since the second term in (|2.16p involves products of local operators at different 
points on the boundary. 

We glossed over this point above, but it is important to distinguish between the 
vacua of the interacting and free theories in the expression (12.91) . The Bunch-Davies 
vacuum that appears as the initial state in (|2.9p is the interacting vacuum, and is 
related to the vacuum of the free theory by 

|BD) oc Texp (-i f drjHiS) |BD) free , (2.17) 

up to an overall factor that was absorbed into the measure in (|2.9p . On the other 
hand, the final state (+| is defined the same way in the interacting theory as in the free 
theory; it is annihilated by e 2lU7T a(k) — a'(— k), where a' and a are the Bunch-Davies 
creation and annihilation operators that appear in the interaction picture fields. 



3 Scattering Amplitudes 

There is an equivalent way of viewing the correlation functions of 0± in terms of a 
boundary scattering matrix, along the lines of |16|. I17j. To see this, let us define a pair 
of solutions to the wave equation by the formulas 

cf>i(r,,k) = K ± (r,,k)e- i *- s , (3.1) 

where the bulk-to-boundary propagators are given by 

K±( V ,k) = (±|0±(fc)*(77,-fc)|O) 
2 Tv 7T ri d ' 2 k ±u 



H ±I/ { V k). (3.2) 



Here |— ) is the state that is annihilated by C+. The solutions 03. and c^L start out 
as plane waves that behave as e ikv , and are fast-falling and slow- falling respectively at 
up to a delta function singularity at the point x, 

4 to k) ~ ((-^ + {lk ^±,) ] { ~ r])A± ) (3 ' 3) 



2 Actually, 15. suggests that the marginal operator in question is O+O-, rather than (a^) 2 . 
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Note that our definitions of 4>% and § x _ are different from those in [16]; in that work, 
the bulk-to-boundary propagators were not purely negative frequency at early times. 

Next, as in [10J, we decompose the field operator <!> as = + where 
^±(Vix) ~ (— t?) A± near X + . It is then straightforward to check that 0± can be 
represented (IS i\ creation operator, 

0±(x) = (a|)t = - f dZ» 0|^$±. (3.4) 

Here we are assuming that asymptotes to a free field at I + ; for simplicity, we have 
left out a wavefunction renormalization for <J>. Thus the matrix elements of 0± between 
|BD) and an arbitrary final state (/| are equal to the scattering amplitudes 

(/ino + (^n°-^)i BD )=^n(4) t n( a *) t i BD )- ^ 

i j i j 

For a generic final state (/|, one can not directly apply an LSZ-like argument to 
(|3.5p . although this can be done by modifying the definition of the bulk-to-boundary 
propagators. However, in the special case when the expectation values are of the form 
(±| Y\i 0±(xi)\BD) , we may generalize the standard flat-space techniques in order to 
derive a diagrammatic expansion for the amplitudes, which is analogous to the Witten 
diagram expansion in AdS |19j . To do this, let us first rewrite the amplitude as 

(±lIJO±(^)|BD) = (±|Jji f dS"4^*|BD>, (3.6) 

where we have discarded terms that arise from the commutator [<3? + (t/ c , x), &-(r/ c , y)], 
and therefore do not contribute to fully connected scattering processes. 

Integrating by parts, using the equations of motion, and discarding all surface terms, 
the amplitude then becomes 

(±|JJO±(^)|BD) = J ]Jd d+1 x\V^K±(4,Xi)(MTl[Hx'i)\VV)A^ (3-7) 

i i i 

where the subscript "A" indicates that the external legs have been amputated. This 
provides a simple relation between truncated bulk correlation functions and matrix 
elements of the operators in the CFT. Note that the oscillatory surface terms at early 
times in (|3.7f) have been killed by a contour rotation (or an infrared regulator, from 
the point of view of the CFT). Indeed, K± is purely negative frequency with respect 
to the Bunch-Davies vacuum at the horizon, and only the negative frequency part of 
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the <E> fields survives in the matrix element. Thus the entire integrand behaves as e lkv 
in the infinite past, and vanishes after taking rj i— > rj + ie. 

We have been referring to (13.51) as the S-matrix, but this is a slight misnomer - a 
more accurate term would be the S-vector [18] , since the initial state is fixed to be the 
Bunch-Davies vacuum. It would be interesting to understand how modifications of the 
initial state are reflected in the CFT. 



4 In-In Expectation Values 



So far, we have only considered transition amplitudes from the Bunch-Davies vacuum 
to some final state. In fact, these amplitudes contain the same information as the 
independent behaviors of Euclidean in-in expectation values near I + (see [20] for a 
comprehensive review of the in-in formalism). To show this, let jP[0+] and G[0_] be 
functions of 0±. Inserting a complete set of states and using (|2.14p . we may express 
the expectation value of the product of F and G as 



(4.1) 



(F[0 + ]G[0-]) BD = I [#+] <BD|^ + )(0 + |F[O + ]G[O_]|BD 
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[d<t> + ]Z* CFT [<t> + \F 



(G[2iv<J )+ ]ZcFT&+}) 



where Zcft [</>+] is the CFT partition function with the source </>_ turned off. By (|2,5p . 
the expectation values (|4.ip determine the behavior of expectation values of <& near the 
boundary, so this behavior is in turn completely determined by the amplitudes (|2.14p . 
Note that insertions of C_ and + respectively encode the leading and subleading 
independent behaviors of expectation values at late times. 

The relation (|4,ip does not treat 0+ and C_ on an equal footing; expectation values 
of C + are computed by functionally differentiating, while expectation values of 0_ are 
evaluated by inserting a 0+ field. To remedy this, let us Legendre transform the CFT 
partition function, 



Zcft[ 



/> 



exp 



-2iv 



d d k 
(2tt)° 



ci>+{k)<j)-{-k) Z CFT [ 



(4.2) 



Plugging this into (14. ip . we then find the simple relation 



(F[0 + ]G[0.]) BB 



/> 



F 



G 



5<j) 

X Z£ FT [(j>+\ Z C FT[(f> 



exp — 2iv 
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which is more convenient for computations than (|4.ip . So we see that Bunch-Davies 
expectation values of 0± are computed by a system of two fluctuating fields <j>±. A 
related approach to computing extrapolated Bunch-Davies expectation values is the 
doubled boundary theory of [21]. 

Note that the Legendre transform of the CFT partition function has a natural 
interpretation. In the AdS/CFT context, Legendre transforming the partition function 
is accomplished by switching to alternate boundary conditions on the bulk field [22J. 
Analytically continuing the resulting bulk correlation functions to de Sitter space gives 
the transition amplitudes 



r <j>( v =o)~(- v ) A - 

{-\T$(x 1 )---$(x n )\BD) = / [#]<Kx!) exp(iS[</>]). (4.4) 

J cf>(r)=~oo)^e lkr > 



Thus the Legendre transform Zcft[<^-] is the generating functional for expectation 
values of C_ between (— | and BD). 

In order to see how the above procedure works in practice, let us take a cubic bulk 
interaction, 



S ^ = -yJ d + x V=9P- (4-5) 

The terms in the wavefunction in (j4.3|) that contribute to the tree-level three-point 
functions are then 

ry u. 1 f f d d k i{ik) 2v T{l-u). A . r . |2 

z CFT [M = ^[j W y d Vr ( ,) l<Mfc)l 



iX 
31 



J d d+1 x' y /^]Jd%K + (x\x i )cj> + (x i U , (4.6) 

i=l / 



and the Legendre transform Zqft is identical but with u, K + , and <p + replaced with 
— v, K-, and </>_. Performing the integrals in (|4.3h . we find the CFT correlators 



(0 ± (k)0 ± (k')) BB = - r(1 ^ ff fc V mfc + fc'), (4.7) 
(0+(jfe)0_(jfe > )) B D = v(i - cot(vir))(2Tr) d 5 d (k + k'), 



\i=l I BD 



Re(-a J |0_ ai; (fc i )0_ ai (-^)|BD) 



where each of the dj's can be equal to either + or — . 

On the other hand, we may compute the two-point and three-point functions di- 
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rectly in the in-in formalism, 



($(77, k)$(rf, A:'))bd = W ^ HZ( V k)H u (r/k)(2ir) d 5 d (k + k'), (4.8) 



where Gbd is the Euclidean Wightman function. Pulling these results to the boundary 
then gives the same answers (|4.7p . 

5 Anomalies 

In this final section, we will investigate a phenomenon that was pointed out in [T| I23j: 
the term in the Hartle-Hawking wavefunction that is proportional to the Weyl anomaly 
drops out of the late-time probablity measure. 

To illustrate this point, let us consider three-dimensional gravity with a positive 
cosmological constant. We may compute the wavefunction by analytically continuing 
the partition function from AdS [3J. Since the central charge of the dual CFT to AdS% 
gravity is equal to 3£Ads/(2G) [23], where Iaas is the AdS radius, the AdS partition 



function takes the forrro 

Zb] = exp (-|^| J d 2 x^d 2 y^R{x)U^R{y)) , (5-1) 

where R and □ are the Ricci scalar and Laplacian formed out of the boundary metric 
7. Here we are assuming that 7 is topologically trivial, so that it is conformal to the 
flat metric 5 a b- Performing the analytic continuation ^Ads = ^ds on (15.ip . we find the 
Bunch-Davies wavefunction 

*bd[7] = exp J d*x^jd 2 y^jR(x)n- l R{y)} . (5.2) 



We see from (|5.2p that the conformal anomaly cancels out when computing the late- 
time probability measure I^bdMI 2 - This implies that if matter fields are added to 
the theory in order to obtain a dynamical metric, then bulk expectation values are 
protected from the conformal anomaly. 

In fact, this observation extends to the matter anomaly as well. For instance, the 
conformal anomaly in AdS/CFT for a bulk scalar corresponding to an operator of 



3 1 would like to thank D. Harlow for explaining this argument. 
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dimension A + = d/2 + n for n a positive integer is [27] 

A[<f>] oc 4dk" / ^jA; 2 «^(z c ,fc)| 2 , (5.3) 

up to a real numerical constant. Under the analytic continuation (z c , £\ds) ^ ( — ^Vc ^ds)> 
this term picks up an overall imaginary factor of i(—l) d ~ n+1 , implying that the anomaly 
again cancels out after squaring the wavefunction. 

Recall that the bulk interpretation of the conformal anomaly is that infrared di- 
vergences break the invariance of the bulk theory under the de Sitter isometry group 
SO(d + 1,1). These divergences are consequences of the fact that de Sitter space has 
infinite volume, and should not be confused with the well-known IR divergences for a 
massless scalar in de Sitter space, which arise from the integration over zero modes in 
the propagator [28]. Thus, at least in the cases that we have considered here, late-time 
Euclidean expectation values of bulk fields are protected from the infinite-volume IR 
divergences in the wavefunction that break the de Sitter symmetries. This restoration 
of de Sitter invariance has been confirmed in a simple model of inflation |29] . 

It is natural to ask if expectation values of the conjugate momentum to (j) are 
also protected from the conformal anomaly. To see that this is not the case, note 
that late-time expectation values of the conjugate momentum II can be computed by 
differentiating the wavefunction with respect to (f>o and integrating over </>q, 

(F[n(r ?c ,x)]) BD = J [d<p ] |*BD[</>o,r? c ]| 2 F fe^^]^^^) ■ ( 5 - 4 ) 

This follows directly from Hamilton-Jacobi theory; the variation of the on-shell action 
with respect to the final boundary condition is equal to the conjugate momentum 
evaluated at the boundary. It is clear from (|5.4j) that the phase of the wavefunction does 
not drop out of expectation values of II. Equivalently, the phase of the wavefunction 
contributes to Bunch-Davies expectation values of <f> fields at different times, as noted in 
PQ. Naively this implies that a metaobserver |30j can compute the conformal anomaly 
by measuring late-time expectation values of II. It would be interesting to understand 
if this effect is indeed physical, or if it can be eliminated via field redefinitions. 

6 Conclusion 

We have shown that computations in perturbative quantum field theory in de Sitter 
space can be reorganized as calculations in a hypothetical conformal field theory living 
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at X + . One immediate application of our results is at a practical level: many of the 
difficult integrals involved in computing AdS/CFT correlators have been evaluated in 
the literature, and these can be carried over to Euclidean expectation values in de 
Sitter space by analytically continuing and using the formalism of Section UJ In addi- 
tion, we expect that the formula (|4.3p is helpful for constraining possible shapes of the 
independent behaviors of late-time expectation values of massive fields, since the con- 
straints imposed by conformal invariance are most easily analyzed from a wavefunction 
perspective [51] . 

Another avenue for future investigation would be to better understand the relation 
between our approach and that of [9] [10l [TTj . In those works, the bulk theory is dual to 
two different field theories on the boundary, depending on whether one imposes Dirich- 
let or Neumann boundary conditions on the bulk field at X + . Here we have chosen 
to impose future boundary conditions by fixing the final state rather than by fixing 
the asymptotic behavior of the bulk field, since this makes the connection to standard 
field theory in de Sitter space more transparent; however, it would be interesting to 
understand the latter perspective as well. 

On a broader note, any complete treatment of de Sitter holography must address 
the question of the finiteness of the entropy (see [52 J for a review) . The bulk dual of our 
CFT Hilbert space is infinite-dimensional in perturbation theory, since one can insert 
any number of operators at arbitrary locations on the boundary. Recent investigations 
suggest that nonperturbative effects may reduce the number of states in the bulk [10J , 
but it is clear that much work still remains to be done. 
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